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Abstract 

We study, by means of three-dimensional Monte Carlo simulations, the current- 
voltage (IV) characteristics and the voltage noise spectrum at low temperatures 
of driven magnetic flux lines interacting with randomly placed point or columnar 
defects, as well as with periodically arranged linear pinning centers. Near the depin- 
ning current J c , the voltage noise spectrum S(uj) universally follows a l/oj a power 
law. For currents J > J c , distinct peaks appear in S(u>) which are considerably 
more pronounced for extended as compared to point defects, and reflect the spatial 
distribution of the correlated pinning centers. 
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The dynamics of flux lines in the mixed state of high-T c superconductors has 
been a topic of intense interest in the past few years [1] . The area of research 
is motivated both by important technological questions and a quest of un- 
derstanding complex non-equilibrium behavior of driven systems. Defects in 
the superconducting materials can either be localized (point pinning centers, 
such as oxygen vacancies) or extended. For example, artificial columnar pins 
are introduced via high-energy ion irradiation to pin the magnetic vortices in 
order to reduce energy loss due to dissipation when current is passed through 
the sample [1]. The non-equilibrium dynamics of flux lines in the presence 
of defects generates non-trivial steady states which can be probed in experi- 
ments measuring, e.g., the current- voltage (IV) characteristics [1], the voltage 
noise power spectrum [2-4], or local flux density noise [3]. Analytical studies 
of the vortex dynamics are usually limited to asymptotic regimes [1,5,6] when 
the current J is either much smaller or much larger than the critical cur- 
rent J c . The intermediate current regime has been investigated by numerical 
techniques such as Langevin molecular dynamics simulations [7-10] or Monte 
Carlo techniques [11,12]. These numerical simulations are largely restricted 
to two dimensions [7-10], while very few three-dimensional studies have been 
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reported [11-13]. Here we present a fully three-dimensional Monte Carlo simu- 
lation of vortex dynamics in the presence of both columnar and point defects. 
We focus on the effects of defect correlations in the voltage noise spectrum 
and the IV characteristics, some of which are inaccessible by two-dimensional 
simulations. 



We study the simplest possible model where the flux lines are non-interacting. 
This may be viewed as a starting point towards a more realistic model that 
takes their mutual repulsion into account. At any rate, our present study 
should apply to a highly dilute vortex system. In spite of the simplistic non- 
interacting flux line model we find non-trivial signatures of broad-band noise 
(BBN) [4,2,3] at J « J c as well as pronounced peaks in the voltage noise 
spectrum in the presence of extended disorder when J > J c . We would like 
to emphasize that our results may apply quite generally for the dynamics 
of extended objects, such as charge density waves [14,15], interfaces [14-16], 
and polymers [14,15], driven through a disordered medium. Hence we con- 
sider a single magnetic vortex in the London approximation as an elastic line 
with line tension energy El[t(s)} = | Jq ds [T~ 2 (dr±(s)/ds) 2 + (dz(s) / ds) 2 }, 
where r(s) = [r_i_(s), z(s)] describes the configuration of the flux line in three 
dimensions, and the crystalline axis c of the superconducting material (as 
well as the external magnetic field) is oriented parallel to z (in a sample 
of thickness L). This linear elastic form of energy holds good as long as 
\dr±(z)/dz\ < 1/r [17], where T denotes the anisotropy ratio [1,5]. The line 
stiffness is given by e « e ln(A ab /£ a6 ) [1,5], with e = (0 o /4vrA afe ) 2 (0 O = hc/2e 
is the magnetic flux quantum). \ ab and (, a b denote the penetration depth 
and the superconducting coherence length respectively in the ab plane. For 
high-T c materials, r > 1 [1], giving the quadratic form of £x,[r(s)] a wider 
range of validity. In the presence of (point or columnar) defects the flux line 
is described by the free energy Fl(T) defined through exp[— /3Fl(T)L] = 
/£>[r(s)]exp [-(3E L [r(s)\ - (3 J L dsV p [r(s)]\ , where (3 = (k B T)-K We model 
the pinning potential as a sum of N p independent potential wells, V^,[r(s)] = 
J2k=i UQ{r p — \r(s) — rj^|), where r p is the radius, &(x) denotes the Heaviside 
step function, and the rjf^ indicate the location of the pinning sites. 



In our Monte Carlo (MC) simulation each flux line is modeled by N = L/a^ 
beads where the ith bead is located at = (xi,yi,Zi) and interacts with its 
nearest neighbors via a simple harmonic potential | J2(j)=i-i, »+i [r~ 2 ( r ±« — 
r ±{j)) 2 + (\ z i ~~ z {j) \ ~ °o) 2 ]- The line is placed in a box of size L x x L y x L z 
with periodic boundary conditions in all directions. The columnar pins and 
spherical pins are respectively modeled by constant cylindrical or spherical 
potential wells of (uniform) strength U and radius b . The cylindrical wells are 
oriented parallel to the c axis. We investigate the dynamics for three different 
distributions of defects, namely, (i) columnar pins distributed randomly or (ii) 
arranged in a square lattice in the xy plane, and (iii) point pins distributed 
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uniform randomly in the sample. When an external current J = — Jy is applied 
it produces a Lorentz force fj, = (l/c)£ x J = (J/c)x, per unit length of 
the flux line. Since in the presence of weak currents J the flux line locally 
moves with an equilibrium dynamics one may incorporate the effect of the 
force in the MC by introducing an additional work term — fj, • J L dsr(s) in 
the energy E L [r(s)} [5,11,12]. At t — the flux line starts off from a straight 
line configuration parallel to the z axis. We have checked that the results in 
the steady state are independent of the initial configurations. At each trial 
a randomly chosen point on the line is updated according to a Metropolis 
algorithm [18]. The point can then move in a random direction a maximum 
distance of A < bo/^/3 to guarantee interaction with every defect. In the 
simulation the values Zi are held fixed; we have checked [19] that results do 
not change qualitatively even if these positions Z{ are allowed to fluctuate. 

The drift velocity of the flux line is proportional to the average velocity of the 
center of mass (CM) v cm = ([R cm (r) — R cm (0)]/r), where R 

emir} -R-cm(O) is 

the distance traversed by the CM in a time interval of r; (• • •) and the overbar 
denote the average over Monte Carlo steps (MCS) in the steady state and over 
different realizations of disorder, respectively. The voltage drop measured in 
experiments is caused by the induced electric field E = B x v cm /c [20]. All 
length and energy scales are measured in units of bo and eo, respectively. The 
average distance between the pins for a uniform random distribution or the lat- 
tice constant for the periodic array of columnar pins is taken as d = 15&o- The 
parameters \ a b, (, a b, £, U, and T are chosen to be 566 , 0.646 , 4e , 0.0075e , 
and 16 respectively in the simulation; these numbers are consistent with the 
experimental data for high-T c materials [1]. Since here we are interested in 
studying the effects of defect correlations in the dynamics, we restrict our 
simulation to temperatures T/T* < 1, where T* = \J eUb§ is the temper- 
ature above which entropic corrections due to thermal fluctuations become 
relevant for pinned flux lines [5]. Therefore thermally induced bending and 
wandering of the flux lines are largely suppressed, and we can interpret our 
results in terms of low-temperature dynamics. In the simulation all the data 
have been collected in the the steady state where the system arrives at t > 10 5 
MCS and at t > 10 6 MCS for columnar and point pins, respectively. The size 
of the system ranges from 60 x 9 x 60 to 2000 x 9 x 60 in the simulations, 
and the data are averaged over 20 to 50 realizations of disorder. The value of 
t ranges from 30 to 500 MCS in the simulation. 

The IV characteristics displayed in Fig. 1 were computed for random distri- 
butions of point and columnar pins, as well as for columnar defects arranged 
in a square lattice in the xy plane at a temperature T = 0.25 x 10~ 3 << T*. 
The critical current J c , which at zero temperature would indicate the loca- 
tion of a non-equilibrium depinning transition, is higher for the columnar 
pins compared point disorder because of the increased defect correlations. 
The value of J c is the same for both random and periodic extended defect 
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Fig. 1. IV characteristics for point pins (*) and columnar defects, both with ran- 
dom (o) and periodic (A) spatial distribution. Note that J c is almost two times 
larger for columnar defects compared to point pins. For currents J > J c the ran- 
domly distributed columnar pins show a larger resistance than those with a periodic 
distribution (see text). 

distributions. The flux line remains localized at few defect sites for J < J c 
as thermal creep happens extremely rarely at T << T* in the simulations. 
For currents J > J c , the line becomes depinned and we observe that for 
columnar defects the ensuing voltage is lower for the periodic distribution as 
compared to the random distribution. This phenomenon can be understood 
in the following way. Suppose the energy of a flux line inside a columnar pin 
(say pin 1) of length L is U\ and the energy of the flux line in the neigh- 
boring pin (pin 2, at a distance a from pin 1) where it arrives after a time 
At(12) is U2, then U2 — U\ ~ —Lf^a, which is the effective energy barrier 
between 1 and 2 [5]. The transit time is At(12) oc exp(— /3L/^a), and hence 
the mean velocity v oc aexp(/3L/ia). Let us take a configuration where the 
three pins A, B and C are arranged in a straight line. The distance between 
A and B is 6ab = a + 5 and the distance between B and C is dsc = a — 5. 
Note 5^0 and 5 = respectively for random and periodic distribution of 
the pins. Here the average velocity of hopping between neighboring pins is 
oc \[d AB /kt{AB) +d BC /At(BC)] oc v[cosh(/3Lf L 5) + 5smh([3Lf L 5)} > v. 
Therefore the average velocity of the flux line is less for a periodic defect ar- 
rangement as compared to a random distribution. For periodic arrangement 
of the pins the IV characteristic depends strongly on the orientation of J with 
respect to the lattice direction, for this determines the effective density of 
defects encountered by the moving vortex [19]. 

The velocity (voltage) noise is calculated by computing the velocity fluctua- 
tions about v cm . We evaluate the power spectrum S{uj) = v x {uj)v*{uj), where 
v x {uj) is the Fourier transform of the velocity fluctuation v x = v x — (v x ), with 
v x (t) = \X cm (t + T) — X cm (r)]/r. The data are averaged over 50 realizations of 
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Fig. 2. (a) Plot of the power spectrum S(co) for J = for columnar pins at tem- 
peratures T/T* = 6.5 x 1(T 4 (bottom), 6.5 x 1(T 2 (middle) and T/T* = 3.3 x 10 _1 
(top). At T/T* = 6.5 x 10~ 4 , r g << 6 and the voltage noise is essentially thermal 
(white). At higher temperatures (for T/T* = 6.5 x 1(T 2 and T/T* = 3.3 x 1CT 1 ) 
the flux line starts to 'feel' the potential well, rendering the fluctuations over large 
length scales correlated, as reflected in the nonlinearity in the small u> behavior of 
S(oj). The nonlinearity follows a uj 2 (solid line) behavior for small u> (see text), 
(b) Log-log plot of S(u) vs. frequency uj for both columnar (A) and point(o) defects 
near the depinning threshold at J J c . The solid line, given as a guide, has slope 
-2. 

disorder for random distributions of defects. At T = and for J < J c , the flux 
line is bound to one or more pins depending on the type of defects present in 
the system, and v cm = 0. For temperatures close to T = the vortex is still 
trapped in the defect potential wells, and r 2 = Y,iLi( r ±i~^±) 2 « &o, where 
ri = J2iLi r ±i- Therefore the line hardly tunnels out of the pins and S(u) 
yields a white (thermal) noise spectrum (Fig. 2 (a)). At higher temperatures, 
when r g b , long length scale fluctuations can tunnel out of the defects 
and produce a nonlinear rise as oc uj 2 at small frequencies in S{ui) (Fig. 2(a)). 
This is because S(oo) is an even function, and S(uo = 0) xs (the equality is 
valid at T = and in the L — > oo limit). The non- white nonlinear part of 
S(u) persists upto larger frequencies for columnar pins as compared to the 
point defects because of increased correlations induced in the motion of the 
flux line by the correlated disorder [19]. At temperatures T > T*, r g >> b 
and the velocity fluctuations are dominated by thermal fluctuations; thus the 
high-frequency white-noise part of S(u>) increasingly dominates the spectrum 
(Fig. 2(a)). 

At J J c , S(u) shows al/w" behavior (Fig. 2(b)), where a a M F = 2 
for almost a decade for both point and columnar defects with deviations at 
small and large frequencies. This is to be interpreted as a remnant of the 
zero-temperature depinning transition, representing a non-equilibrium critical 
phenomenon [21-23]. Its analysis through a functional renormalization group 
calculation gives a — 1.5 for point defects to one- loop order in three dimen- 
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Fig. 3. (a) Plot of S(lo) vs. Lo/v cm for columnar pins arranged in a periodic array 
for J = 0.05(o), J = 0.07 (square), and J = 0.09(A). The peaks occur at the 
same scaled frequency u> oc v cm /d. Note that the width of the peak decreases as 
J is increased indicating the wider separation between the time scales r r and t v . 

(b) Plot of S(oj) vs. co for columnar pins distributed randomly in the sample for 
J = 0.05 (square). Note the presence of secondary peaks in the spectrum (see text). 

(c) Plot of S(u>) vs. oj for point pins at J = 0.05 (square). The pronounced peaks 
present in the case of columnar pins are considerably suppressed (see text). 



sions [22,23], with corresponding mean-field value olmf — 2. Thus, while our 
data for lowest frequencies are affected by our finite system size, the uj~ 2 decay 
represents the finite-temperature mean-field like behavior to be expected off 
criticality. The effective exponent a then crosses over to smaller values, per- 
haps even its true critical value, at larger frequencies. In experiments [4], the 
broad-band noise (BBN) features show a u~ a decay with a ~ 2 which may 
similarly be a manifestation of the single-line depinning transition at a finite 
temperature. 

For currents J > J c , two distinct time scales appear in the dynamics — the 
residence time r r which is the duration of the vortices being trapped on the 
defects, and t v , the time taken by the flux line to travel between the pins. 
In the regime where these time scales are well separated and r r < r v , we can 
expect t v to produce pronounced peaks in the velocity noise spectrum. In Fig. 
3(a) we plot S(uj) at different driving currents for the periodic distribution 
of columnar pins which reveals there is a single peak corresponding to the 
frequency Co oc v cm /d oc l/r v . Note that the width of the peak decreases as 
the external current J increases, which widens the separation between r r and 
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t v . In a random defect distribution (Fig. 3(b)) there would be a distribution 
of the time scale r v around the average value d/v cm . Therefore there could be 
many secondary peaks present in the spectrum arising due to the distribution 
of t v in addition to the principle peak which corresponds to r r oc d/v cm . 
The secondary peaks would eventually diminish as the distribution of the 
time scales will narrow down at d/v cm with more number of averaging on the 
defect configurations. Unlike for correlated defects, the time scales r r and r v 
are not well separated for point pins. The plot in Fig. 3(c) reveals, as expected, 
that the peaks are considerably suppressed as compared to the situation for 
columnar defects (Fig. 3(b)). This difference in S(u) between columnar and 
point pins which results due to the lack of spatial correlation in the z direction 
for the point disorder can be observed only in a three-dimensional dynamical 
simulation. It is due to a generic difference between correlated and localized 
pinning potentials, and should survive even in the presence of interactions 
between the flux lines. The ensuing absence or presence of the narrow-band 
noise peaks in the voltage noise spectrum S(uj), and its specific features, can be 
utilized as a signature to identify and characterize the type of disorder present 
and responsible for flux pinning in experimental superconducting samples. 

In conclusion we have investigated the three-dimensional dynamics of non- 
interacting flux lines in the presence of both point and correlated defects. 
We find the nature of the current-voltage characteristics to depend on the 
defect correlation and spatial distribution of the pinning centers. The voltage 
noise spectrum shows pronounced peaks in the case of correlated pins. The 
difference in the power spectrum at J > J c between columnar and point 
disorder can serve as a novel diagnostic tool to characterize the pinning centers 
in superconducting samples. 
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